Universal fluctuations in KPZ growth on one-dimensional flat substrates 
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We present a numerical study of the evolution of heiglit distributions (HDs) obtained in interface 
growth models belonging to the Kardar-Parisi-Zhang (KPZ) universality class. The growth is done 
on an initially flat substrate. The HDs obtained for all investigated models are very well fitted by 
the theoretically predicted Gaussian Orthogonal Ensemble (GOE) distribution. The first cumulant 
has a shift that vanishes as t~^^^ , while the cumulants of order 2 < n < 4 converge to GOE as t~'^^^ 
or faster, behaviors previously observed in other KPZ systems. These results yield a new evidence 
for the universality of the GOE distribution in KPZ growth on fiat substrates. Finally, we further 
show that the surfaces are described by the Airyi process. 
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The last decade witnessed a great advance in the 
understanding of paradigmatic nonequilibrium models 
for interface growth by means of rigorous analytical re- 
sults [nil] as well as their experimental realizations [3111]. 
Among the most prominent examples are several models 
belonging to the Kardar-Parisi-Zhang (KPZ) universality 
class introduced by the stochastic equation [S] 

f = -V2^+^(V/.)^+^, (1) 

that describes the evolution of an interface h{x,t). In 
this equation, is a white noise defined by (77) = and 
{r]{x, t)Tj{x', t')) = DS{x - x')S{t -t'). In d = 1 + 1, the 
KPZ equation yields self-afSne surfaces with an interface 
width in a scale of length £, defined as w = A/(7i^y^r'(/jy2"^ 
given by the Family- Vicsek scaling ansatz [6]: w{£,t) = 
t^J'(€/i^/^), where T{x) ^ cont. for x > 1 and T{x) ^ 
x^^ for X ^ 1. In = 1 + 1, the exponents of the KPZ 
class are /3 = 1/3 and z = 3/2 _5j. Several models with 
the scaling exponents of the KPZ universality class have 
been reported [71[H]. 

For system belonging to the KPZ class in d = 1 + 
1, the Family- Vicsek scaling suggests an interface height 
evolving as pj 

hit)~v^t + sign{X)irty^^X, (2) 

where Wqo £^nd F are non-universal (model dependent) pa- 
rameters, and X is a time independent random variable. 
The height distributions (HDs) and consequently x were 
computed exactly for some models in the KPZ class (see 
Refs. [2 [2] for recent reviews), strongly suggesting that x 
is a universal feature of the KPZ class. In a pioneer work, 
Johansson [S] established a link between random matrix 
theory and KPZ class by determining the exact asymp- 
totic HDs of the totally asymmetric exclusion process 
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(single step model) as the Tracy- Widom (TW) distribu- 
tion of the Gaussian unitary ensemble (GUE) [10] . In the 
same year, Prahofer and Spohn [TTJ [T^] determined the 
HDs of the radial polynuclear growth (PNG) model as a 
GUE distribution whereas the Gaussian orthogonal en- 
semble (GOE) distribution was obtained for the growth 
on a flat substrate. Quite recently, exact solutions of the 
KPZ equation in ci = 1 -I- 1 corroborated GOE distribu- 
tions for flat [13] and GUE for radial [Ml [15] geometries, 
respectively. 

Besides TW distributions for the heights, the limiting 
processes describing the surfaces were identified as the 
so-called Airy processes [1]. For circular growth, ana- 
lytic solutions [16 , experiments [3], and simulations [17] 
strongly suggest that the limiting process ruling the KPZ 
universality class is given by the dynamics of the largest 
eigenvalue in Dyson's Brownian motion of GUE matrices, 
the Airy2 process. Analogously, a limiting process for 
the flat case, named Airyi, was found for PNG [TH] and 
a similar model [I^, both belonging to the KPZ class, 
indicating that it may be a universal feature of the KPZ 
class. It is worth mentioning that Airyi process does 
not correspond to the dynamics of largest eigenvalue in 
Dyson's Brownian motion of GOE matrices [^ . 

KPZ universality class in l-fl dimensions was convinc- 
ingly verified in a meticulous experiment involving two 
turbulent phases in the electroconvection of nematic liq- 
uid crystal films [21 H] ■ Using an unprecedented statis- 
tics for experimental systems (~ 10'^ realizations), it was 
possible to observed not only the KPZ scaling exponents 
/3 = 1/3 and z = 3/2, but also the GUE and GOE dis- 
tributions for radial and flat geometries, respectively. In 
both cases, cumulants of order n = 2, 3 and 4 converge 
fast to the corresponding GOE and GUE cumulants, but 
the mean (n = 1) approaches the theoretical value as a 
power law t^^^^. This decay was also found in an ana- 
lytical solution of the KPZ equation in d = 1 -I- 1 with 
an edge initial condition |14| and in simulations of ra- 
dial Eden model [T7]. However, Ferrari and Frings [5T] 
have shown the non-universality of the amplitudes of this 
correction and that higher order cumulants have no cor- 
rection up to order 0{t~^^^). 



model 




r 


A 


BD 


2.1398(7) 


4.90(5) 


2.7(2) 


RSOS 


0.41903(3) 


0.252(1) 


0.81(3) 


KPZ I 


0.05725(4) 


0.00081(1) 


0.008(1) 


KPZ II 


0.08130(5) 


0.00493(2) 


0.035(2) 



TABLE I: Non-universal parameters for the investigated mod- 
els. 



The analysis of height distributions in computer sim- 
ulations of non-solvable models, supposedly belonging to 
the KPZ class, are relatively less frequent than the an- 
alytical counterpart of solvable models. Indeed, some 
reports for flat geometries dating from the beginning 
of nineties, when computer resources were quite lim- 
ited, demonstrate a surprisingly good agreement with the 
lately found out theoretical GOE distribution [22l |23] . 
For a radial geometry, we have recently observed asymp- 
totic GUE distributions in computer simulations of dis- 
tinct Eden models grown from a single seed [I7j . 

In the present work, we revisit the computer simula- 
tions of non-solvable models having the KPZ exponents 
in 1 -|- 1 dimensions. A flat initial condition is used. 
We also integrated the KPZ equation numerically. In all 
cases, we confirmed that HDs and cumulants agree with 
the GOE distribution. As recently observed for other 
systems [3J HI [T3], a correction in the mean approach- 
ing zero as t~^^^ was found while corrections in higher 
order cumulants consistent with t^'^/^ were obtained for 
three of four investigated models. These results are in 
agreement with the corrections obtained analytically for 
solvable KPZ models [2l]. Finally, the covariance of all 
investigated models agrees with the Airyi process, as ver- 
ified for some solvable models [THl HH] ■ 

We study the restricted solid-on-sohd (RSOS) and bal- 
listic deposition (BD) models as well as the numerical 
integration of Eq. ([T]) for two sets of parameters. In all 
cases, we use one-dimensional lattices of size L = 2^° and 
periodic boundary conditions. Statistical averages were 
performed over at least 400 independent samples. The 
RSOS model consists in the deposition of a particle in 
a randomly chosen site. The deposition is rejected ev- 
erytime the height difference between nearest neighbors 
is greater than one lattice unity. In the BD model, the 
deposition site is randomly chosen and the particle falls 
normally to the substrate. When the first contact with 
the deposit happens, it irreversibly sticks in this position. 

We investigated a discrete one-dimensional KPZ equa- 
tion given by [53] 



+ -(/ij_i - hj+if^ + crVl2Ati?, (3) 
where a = yjQJjfKx and i? is a random number uni- 
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FIG. 1: (Color online) Growth velocity v against t~'^^^ for 
the BD model. Inset shows F against time for BD. Dashed 
line represents the estimate shown in Table |l] 



formly distributed in the interval [—0.5,0.5]. Integration 
instabilities are prevented by replacing the gradient term 
in KPZ equation by f{x) = (1 - e^CI^''l')/C where C 
is a parameter to be adjusted to control the instabil- 
ity . Integration was performed for fixed parameters 
Aa; = 1, At = 0.04 and v = 0.5. Two sets of parameters 
were studied. In the first one, named as KPZ I, we used 
a — 0.1, the coupling constant g = X^D/v^ — 24, and 
C = 0.5. In the second one, the KPZ II, we used a ~ 0.2, 
g = 12 and C = 0.2. 

Comparison between HDs and random matrix distri- 
butions needs highly accurate estimates of the param- 
eters Voc and F given by Eq. In BD and RSOS 
models studied here, estimates are reported in litera- 
ture [22] , but the accuracy is not sufficient for the aims of 
the present work. Thus, we determined the parameters 
for all studied models. Using Eq. ([2]), the asymptotic 
growth velocity can be obtained from the extrapolation 
of plots V = d{h)/dt against t~'^/'^ when t — >• oo, as il- 
lustrated in Fig. [T] for the BD model. The asymptotic 
velocities for all investigated models are presented in Ta- 
bleU Notice that our estimates for RSOS and BD models 
are consistent with those formerly reported in Ref. [22j : 
Voo = 0.419 and 2.14. The parameter F was obtained 
from the second order cumulant of the HDs. Accord- 
ingly Eq. ([2|, it reads as (ft.^)^ ~ {TtY/^{x^)c, where 
(x")c is the nth cumulant of x. Since we expect that x 
fiuctuations are given by the GOE distribution, we used 
{x^)c — {Xgoe)c = 0.63805 [11] in our calculations. The 
insertion to Fig. [l] shows F ~ (('i^)c/(Xgob)c)^^^/^ as 
a function of time for the BD model. The estimated F 
values for all investigated models are presented in Table 
|I] Plots similar to Fig. [Tjwere obtained for all models. 

We study the rescaled HDs following the KPZ ansatz 
given by Eq. If this equation holds, the fiuctuations 
of 9 = (/i - Voot)/{Tty/^ would be given by the GOE 
distribution. Fig. [2]shows the rescaled distributions P{q) 
for all investigated models. As expected, the agreement 
with GOE is noticeable, but the distributions are slightly 



FIG. 2: (Color online) Height distributions scaled accordingly 
Eq. ^ for all investigated models. Simulation time of f = 2 x 
10'* was used for BD and RSOS and t = 800 for the integrated 
KPZ equation. The solid line is the GOE distribution. The 
inset shows a zoom around the peak of the distributions. 
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FIG. 3: Absolute differences between cumulants of the scaled 
HDs and GOE distribution against time. Results for (a) first, 
(b) second, (c) third, and (d) fourth order cumulants are 
shown. Results for KPZ I are shown apart in (e). The dashed 
lines have slopes -1/3 in (a) and (e) and -2/3 in (b)-(d). 



shifted in relation to GOE, as highlighted in the inset of 
Fig.i 

The absolute differences between P{q) and GOE cu- 
mulants, \Kn\, are shown as functions of time in Fig. [s] 
The first cumulant difference decays closely to a power 
law for all investigated models, except for KPZ I. 
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FIG. 4: (Color online) (a) Rescaled local roughness for dif- 
ferent models. Simulation times are f = 5 X 10" for BD and 
RSOS and t = 2 X 10^ for KPZ I and II. (b) Scaled two-point 
correlation function against u = {Al/2)/ (Kt)'^^^ for different 
models. Time is t = 10'' for BD and RSOS and t = 400 for 
the integrated KPZ equation. The solid line is the covariance 
of the Airyi process. 



In the latter, the initial decay is faster but seems to con- 
verge to t^^/^. Ferrari et al. [2l] have shown that the 
corrections in the first cumulant for a number of solvable 
models belonging to the KPZ class read as 



-2/3N 



(4) 



where a is a model dependent constant that can even 
be null for some models. Numerically, if the factor a is 
much smaller than the factor of higher order correction, 
one must observe a fast initial decay [©(t"^/^)] followed 
by a t~^^^ decay. This is observed in the KPZ I results 
shown in Fig. [sj^e). For KPZ I, the higher order cumu- 
lants reach GOE very fast and then the differences |k„| 
fluctuate around zero. Furthermore, the non-universality 
of the constant a is verified in our simulation since a is 
negative for BD and RSOS and positive for the KPZ 
equation integration. Ferrari et al |21| have also shown 
that the higher order cumulants have no correction up to 
order 0{t~'^^^). Our simulations are in agreement with 
this claim as one can see in Figs. |3jb)-[3][e), where some 
curves are consistent with a decay while others de- 

cay faster, which is particularly evident in KPZ I case. 
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The conclusion of Ref. [5T] stating the non-universality 
of the correction refers only to the amplitudes since the 
scenario suggested by Eq. Q is quite general: any cor- 
rection to the first moment, if it exists, should decay with 
an exponent —1/3, while corrections in higher order cu- 
mulants, if they exist, must decay as fast as or faster than 

^2/3. 

Besides the universality of the height fluctuations, the 
limiting process describing the interface is an important 
issue on surface growth. To check if our simulations yield 
surfaces given by an Airyi processes, we studied the two- 
point correlation function 

C{1, t) EE {h{x + I, t)h{x, t)) - (5) 

that is expected to scale as C{l,t) ~ {Kt)^^^gi{u) with 
u = {Al/2)/{Kt)^/^, where gi{u) is the covariance of the 
Airyi process and K = XA^ = 2r (TB]. The amphtude A 
can be obtained from the squared local roughness 

w^l, t) = {h\x, t))i - {h{x, t))f ~ Al/6, (6) 

where (•••); means averages over windows of size I [51l22j. 
Fig. ^a) shows the plots used to determine the ampli- 
tudes shown in Table |lj The short plateau for the BD 
model is due to strong finite time effects of this model. 
Notice that the amplitude is given hy A = D/iy ioi the 
KPZ equation, implying that A = 0.01 and A = 0.04 are 
expected for KPZ I and II, respectively. The slightly dif- 
ferent amplitudes observed are due to the approximated 
integration method . Figure |4]^b) shows the results 



for the rescaled two-point correlation function exhibiting 
a very satisfactory agreement among all models and the 
covariance of the Airyi process (gi). This confirms that 
all models are described by the Airyi process and gives 
one more evidence of this universal feature in KPZ uni- 
versality class. Similar results were obtained at different 
times, indicating small finite-time effects in the two-point 
correlation function. 

In conclusion, we investigated non-solvable models in 
the KPZ universality class in one-dimensional flat sub- 
strates as well as a numerical integration of the KPZ 
equation. We analyzed the height distributions and the 
two-point correlation function. Height fluctuations of 
all investigated models are described by GOE distribu- 
tions as conjectured for the KPZ class for flat substrates, 
and previously observed in analytical jl3J and experimen- 
tal [4] works. The cumulant analysis yields a correction 
vanishing as t~^^^ for the first cumulant whereas the cor- 
rections for higher order cumulants go to zero as i^^/"^ or 
faster. The latter result is in agreement with theoretical 
analysis of Ferrari et al pi] . Finally, we have found that 
the limiting processes describing the surface is the Airyi 
process in agreement with previous results for solvable 
models [18l[T9]. 
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